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Abstract 

Grossglauser and Tse (2001) introduced a mobile random network model where each node moves independently 
on a unit disk according to a stationary uniform distribution and showed that a throughput of 6(1) is achievable. 
El Gamal, Mammen, Prabhakar and Shah (2004) showed that the delay associated with this throughput scales as 
9 (n log n), when each node moves according to an independent random walk. In a later work, Diggavi, Grossglauser 
CI and Tse (2002) considered a random network on a sphere with a restricted mobility model, where each node moves 

(3JT)' along a randomly chosen great circle on the unit sphere. They showed that even with this one-dimensional restriction 

on mobility, constant throughput scaling is achievable. Thus, this particular mobility restriction does not affect the 
' throughput scaling. This raises the question whether this mobility restriction affects the delay scaling. 

This paper studies the delay scaling at 8(1) throughput for a random network with restricted mobility. First, a 
variant of the scheme presented by Diggavi, Grossglauser and Tse (2002) is presented and it is shown to achieve 
8(1) throughput using different (and perhaps simpler) techniques. The exact order of delay scaling for this scheme 
is determined, somewhat surprisingly, to be of 8(nlogn), which is the same as that without the mobility restriction. 
Thus, this particular mobility restriction does not affect either the maximal throughput scaling or the corresponding 
delay scaling of the network. This happens because under this 1-D restriction, each node is in the proximity of 
O . every other node in essentially the same manner as without this restriction. 

Index Terms 

> 

t-j. ' Random wireless networks, scaling laws, constant throughput scaling, delay, 1-D mobility. 

OO ■ I- Introduction 

^ ' Gupta and Kumar [7] introduced a random network model for studying throughput scaling in a fixed wireless 
Q . network (that is, when the nodes do not move). They defined a random network to consist of n nodes where each 
""^5 ' node is distributed uniformly and independently on the unit sphere in M 3 . The network has n/2 distinct source - 
O ■ destination pairs formed at random. Each node can transmit at W bits-per-second provided that the interference 
I> . is sufficiently small. They showed that in such a random network the throughput scales as Q(l/\/n\ogn) per 
^> | source-destination (S-D) pair. 

Grossglauser and Tse [8] showed that by allowing the nodes to move, the throughput scaling changes dramatically. 
Indeed, if node motion is independent across nodes and has a uniform stationary distribution, a constant throughput 
scaling (0(1)) per S-D pair is feasible. This raised the question: what kind of mobility is necessary for achieving 
constant throughput scaling? Diggavi, Grossglauser and Tse [3] considered a restricted mobility model where 
each node is allowed to move along a randomly chosen great circle on the unit sphere with a uniform stationary 
distribution along the great circle. They showed that a constant throughput per S-D pair is feasible even with 
this restricted mobility model. Thus they established that node motion with a stationary distribution on the entire 
network area is not necessary for achieving constant throughput scaling. 

El Gamal, Mammen, Prabhakar and Shah [5] (see [6] for complete details) determined the throughput-delay 
trade-off for both fixed and mobile wireless networks. In particular, it was shown that for mobile networks at 
throughput of 0(1), the delay is 0(nlogn). For mobile networks, the mobility model consisted of each node 
moving independently according to a symmetric random walk on a yjn x yjn grid on the unit torus. 

The constant throughput scaling result of [3] for a network with restricted mobility raises the question whether the 
high throughput in spite of restricted mobility is at the expense of increased delay. Motivated by this question, we 
study the delay scaling for constant throughput scaling in a network with restricted mobility. Somewhat surprisingly, 
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we find that delay scaling is not affected by this mobility restriction either. That is, delay scales as 0(n log n), which 
is the same as the delay scaling when mobility is not restricted. This paper is a consolidation of the preliminary 
work presented in [10]. 

This seemingly surprising result can be explained as follows. Since there are n nodes in a network of constant 
area, the neighborhood of each node is 6(l/n). Based on this, let us say that two nodes meet or are neighbors when 
they are within a distance of The following condition ensures constant throughput scaling in the mobile 

network models presented in [8], [3] and this paper: for 0(l/n) fraction of the time, each node is a neighbor of 
every other node with only ©(1) other nodes in its neighborhood. This ensures that the total network throughput is 
0(n) and that it is distributed evenly among the n/2 S-D pairs, so that the throughput is 0(1). Delay is determined 
by the first and second moments of the inter-meeting time of the nodes. In the case of unrestricted mobility, the 
inter-meeting time of any two nodes is equivalent to the inter-visit time to state (0, 0) for a 2-D random walk on a 
y/n x yfn grid. In the restricted mobility case also the inter-meeting time turns out to be equivalent to the inter-visit 
time to state (0, 0) for a slightly different random walk. However the first two moments are still of the same order 
and hence the queueing delay is the same, leading to the same delay scaling. As a result, even with this particular 
mobility restriction, the maximal throughput scaling and the corresponding delay scaling remain unchanged. 

The rest of the paper is organized as follows. In Section |H] we introduce the random mobile network model, 
some definitions and notation. In Section we present a scheme using random relaying and show that it achieves 
constant throughput scaling. In Section |Wj we show that the delay for this scheme is 0(nlogn) using results 
which are proved in Section |V] The proof of delay of 0(nlogn) consists of analyzing a queue at a relay node in 
two parts. The first part presented in Section llVl identifies an i.i.d. component that is embedded in the arrival and 
service processes of the queue. The second part breaks the dependence between the arrival and departure processes 
by introducing a virtual Bernoulli server. The queueing analysis that follows is carried out in Section |V] 

II. Models and Definitions 

In this section, we present the network model, and the definitions of the performance metrics - throughput and 
delay. We begin by reminding the reader of the order notation: (i) f(n) = 0(g(n)) means that there exists a constant 
c and integer ./V such that f(n) < cg(n) for n > N. (ii) f(n) = o(g(n)) means that lim^oo f{n)/g{n) = 0. (iii) 
f(n) = n(g(n)) means that g(n) = 0(f(n)), (iv) f(n) = u{g(n)) means that g(n) = o(/(n)). (v) f(n) = 0(/(n)) 
means that f(n) = 0(g(n)); g(n) = 0(f(n)). 

Now let us recall what is meant by the uniform distribution of great circles on a sphere. Let S 2 denote the 
surface of a sphere in M 3 with unit area. For x G S 2 , let x' G S 2 be the diametrically opposite point of x. Let 
G{x) denote the great circle obtained by the intersection of S 2 with the plane passing through the center of S 2 and 
perpendicular to the line xx'. Let x be called the pole of G(x). If the pole of a great circle is chosen according to 
a uniform distribution on S 2 then the great circle is said to have a uniform distribution. 

Definition 1 (Natural random walk): A natural random walk on a discrete torus of size m is the process S(t) G 
{0, . . . , m — 1}, t = 0, 1, . . ., such that S(0) is uniformly distributed over {0, . . . , m — 1} and S(t + 1) is equally 
likely to be any element of {S(t), S(t) — 1 mod m, S(t) + 1 mod m}. 

This differs from a simple random walk, where S(t + 1) is equally likely to be any element of {S(t) — 1 
mod m, S(t) + 1 mod m}. Since we are interested only in scaling results, we use the terms (simple) random walk 
and natural random walk interchangeably. 

Definition 2 (Random network): The random network consists of n nodes that are split into n/2 distinct source- 
destination (S-D) pairs at random. Time is slotted for transmission. Associated with each node is a great circle of 
S 2 chosen independently according to a uniform distribution. 

The great circle of each node has ^fn equidistant lattice points numbered from to y/n— 1 placed on it arbitrarily 
resulting in a one-dimensional discrete torus of size yjn. Each node moves according to a natural random walk 
on these lattice points on its great circle. Figure [2 shows a realization of the random network model. Note that 
since the sphere has unit area, its radius is 1/2^. Hence each great circle has perimeter ^pn because of which 
the distance between two adjacent lattice points is yjir/n. 

Let the distance on the sphere between nodes i and j be denoted by d(i,j). We assume the Relaxed Protocol 
model [5] similar to the Protocol model in [7] for successful transmission. 
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Fig. 1. A realization of the random network model. Only the lattice points on the great circles of nodes i and i are shown. The intersection 
of their great circles is Zij. The shaded circle is dj and i and j become neighbors when they are at the two dark lattice points. 



Definition 3 (Relaxed Protocol Model): A transmission from node i to node j is successful if for any other 
simultaneously transmitting node k, 

d(k,j)>(l + A)d(i,j) 

for some A > 0. If a transmission is successful then communication occurs at a constant rate of W bits-per-second. 
For simplicity, we assume that time-slots are of unit length so that when a successful transmission occurs a packet 
of size W is communicated. 

In the other commonly used model (e.g., [8], [7], [3]), known as the Physical model, a transmission is successful 
if the Signal to Interference and Noise Ratio (SINR) is greater than some constant. It is well known [7] that the 
Protocol model is equivalent to the Physical model when each transmitter uses the same power. 

The differences between this model and the model in [3] are: (i) the Relaxed Protocol model is used instead of 
the Physical model, and (ii) each node is assumed to move according to a natural random walk instead of just a 
stationary, ergodic motion with uniform stationary distribution on the great circle. However, this model has the same 
1-D mobility restriction. Further, the proofs clearly show that the assumption of mobility according to a natural 
random walk is not necessary for achieving constant throughput scaling and is used only for computing delay. 

Definition 4 (Scheme): A scheme IT for a random network is a sequence of communication policies, (IT n ), where 
policy IT„ determines how communication occurs in a network of n nodes. 

Definition 5 (Throughput of a scheme): Let Bn n (i,t) be the number of bits of S-D pair i, 1 < i < n/2, 
transferred in t time-slots under policy IT ra . Note that this could be a random quantity for a given realization 
of the network. Scheme n is said to have throughput Th(n) if 3 a sequence of sets An(n) such that 

Au{n) = : min lim inf -Bu (i,t) > Tn(n)l 

[ \<i<n/2 t^co t J 

and P(Aj\(n))^\ as n^oo. 

We allow randomness in policies. Hence, P (Ajj(n)) denotes the probability of Au(n) over the joint probability 
space that captures randomness in the policy as well as the random network instance. We say that event A occurs 
with high probability (whp) if P(A)^l as n^oo. 

Definition 6 (Delay of a scheme): The delay of a packet is the time it takes for the packet to reach its destination 
after it leaves the source. Let D^^j) denote the delay of packet j of S-D pair i under policy U n , then the sample 
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mean of delay for S-D pair i under LI n is 

k 

^=lmsupi££>k(j). 

fc— >oo ™ ■ i 

The average delay over all S-D pairs for a particular realization of the random network is then 

re/2 

Sn„ = -E^„. 

1=1 

The delay for a scheme II is the expectation of the average delay over all S-D pairs, i.e., 

2 n/2 

D n {n)=E[D nn ] = ~Y d E[D\ ln ]. 

i=l 

Now observe that some realizations of the random network may result in the configuration of nodes being such 
that it is not possible to achieve constant throughput scaling. Hence we first define a typical configuration which 
captures the fact that the distribution of great circles is sufficiently uniform everywhere on the sphere. We need 
some notation to introduce this definition. 

Let Gi denote the great circle of node i € {1, . . . , n}. For any two nodes % ^ j, Gi and Gj are not identical with 
probability 1 under the random network model. Two distinct great circles must intersect in exactly two points. For 
each pair % ^ j, select one of the two distinct intersection points of Gi and Gj uniformly at random and call it Zij. 
Let Cij denote the disk on the sphere centered at with radius (2 + A)y/ir/n. See Figure [2 for an illustration. 

Definition 7 (Typical configuration): A configuration (i.e., realization of the random network) is said to be typical 
if the number of great circles passing through each Cij is 6 (\/n). 

Definition 8 (Neighbor): We say that nodes i and j are neighbors at time t if both nodes i and j are at the 
lattice points of their respective great circles that are closest to z%j. 

In Figure [Q the lattice points for nodes % and j that are closest to have been darkened. Under the random walk 
model, it is possible that in some time-slot, a node may not have any neighbors. 

III. Scheme with Constant Throughput Scaling 

In this section we present Scheme II and show that it achieves constant throughput scaling. In the next section its 
delay scaling will be analyzed. Before presenting the scheme, we prove a property of the random network model 
which makes the scheme feasible. 

Lemma 1: Configurations are typical whp. 

Proof: Consider any two nodes i and j. First note that the probability that Gi and Gj coincide is zero. Also 
any two distinct great circles necessarily intersect at exactly two points. By definition, Cij has area c\/n since it 
has radius (2 + A)y7r/n. 

Let ifc, k = 1, ... ,n, k 7^ i,j, be an indicator random variable for the event that the great circle of node k, 
Gk, passes through Cy. By definition, are i.i.d. Bernoulli random variables with parameter p, where p = cij^fn 
where C2 is a positive constant. This is because a great circle passes through a disk of radius it! if and only if 
its pole lies in an equatorial band of width 2R. The probability of this event is Q(R) as the position of pole is 
uniformly distributed over the sphere. 

Thus, the total number of great circles passing through Cij is given by a random variable X = ^ fc 1^ with 
E[X] « O.Sciy^ = 0(y / n). An application of the well-known Chernoff bound for the sum of i.i.d. Bernoulli 
random variable (e.g., see [11]), yields 

P{\X - E[X}\ > 5E[X}} < 2exp (-5 2 E[X}/2) 

= J,, for^V 2(10g g| OSn) - W 

The choice of 8 in © shows that X < c\^fn or X = (\/n) with probability at least 1 — 1/n 3 . Hence by the 
union bound over all n(n — l)/2 possible Cij for i, j = 1, . . . ,n, we obtain that with probability at least 1 — 1/n, 
the number of great circles passing through each Cij is {y/n). ■ 
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A. The Scheme 

The operation of Scheme II depends on whether the configuration is typical or not. If the configuration is not 
typical, direct transmission is used between the S-D pairs along with time-division multiplexing. That is, the sources 
transmit to their destinations once in 2/n time-slots in a round-robin fashion. If the configuration is typical then 
Policy S n as described below is used. Policy S n is a variant of the policies presented in [8], [3]. 

Policy S n : 

1) Each time-slot is divided into two sub-slots - A and B. 

2) Sub-slot A 

(a) Each source node independently becomes active with probability p& > 0. 

(b) If an active node has one or more neighbors then with probability < a < 1, it chooses one at random 
and a packet intended for its destination is transmitted to this randomly chosen neighbor, which acts as 
a relay node. 

3) Sub-slot B 

a) Each node independently becomes active with probability p& > 0. 

b) If an active node has one or more neighbors that are destination nodes, it chooses one at random. The 
active node, which acts as a relay, transmits a packet intended for this destination node, if it has any, 
in FIFO order. 

In policy E n , each node acts as a relay for all the other n/2 — 1 S-D pairs. A packet reaches from its source to its 
destination as shown in Figure |2] A source node, S, transmits its packet to a random relay node, R, which may also 
happen to be the destination itself. The random relay node then moves around carrying the packet. Finally, when it 
becomes a neighbor of the destination, D, the packet is transmitted to D. A relay node may receive several packets 
from a source before it gets a chance to transmit to the destination. To handle this, each relay node maintains a 
separate queue for each of the other n/2 — 1 S-D pairs. 




Fig. 2. Source node, S, transmits its packet to a random relay node, R. The packet is carried by R, until its transmission to the destination 
node, D, when R and D become neighbors. The dotted great circles correspond to other nodes which can act as relays. 

The actual mechanism is slightly more complicated. Since each node decides to transmit at random, it is possible 
that two nearby nodes transmit simultaneously so that transmission is not successful under the Protocol model. In 
order to analyze the throughput of Scheme II, we first state a result about the probability of successful transmission 
between two nodes when they are neighbors under policy S n . 

Lemma 2: Under policy S n , the following hold in a typical configuration. 



6 



(a) In sub-slot A, if nodes S and R are neighbors of each other, S transmits a packet to R successfully with a 
strictly positive probability, independent of n. 

(b) In sub-slot B, if nodes R and D are neighbors of each other, R transmits a packet to D successfully with a 
strictly positive probability, independent of n. 

Proof: We shall only prove for the case of sub-slot A since the proof for the other part is similar. Consider 
a sub-slot A in which S and R are neighbors. Let E\ be the event that S becomes active and E2 be the event that 
S chooses R as a random relay and no other source node in Csr becomes active. If both events E\ and E2 occur, 
S transmits to R and the transmission is successful under the Relaxed Protocol model. Thus, 

P( S transmits to R successfully) = P(Ei n E 2 ) 

= P(E l )P(E 2 \E l ). (2) 

From the description of Policy S n it is clear that P(E\) = apA, which is a strictly positive constant. Next we 
compute P(E<i\Ei) and show that it is lower bounded by a strictly positive constant, independent of n, which will 
imply the statement of the lemma. 

Given that S is active, the probability of successful transmission to R depends on how many other nodes are 
present in Csr since these nodes could interfere, i.e., transmit simultaneously so that the transmission from S to R 
is not successful under the Relaxed Protocol model. 

Since we have a typical configuration, Q(^Jn) distinct great circles or source nodes intersect Csr- Moreover each 
great circle has 6(1) lattice points that are in Csr- For a natural random walk on a discrete torus of size ^fn, the 
probability of being at any particular position is l/yfn. Hence the probability that any of the (\fn) source nodes 
whose great circles intersect Csr is present in Csr with probability ®{l/^/n). Due to the independent movement 
of all nodes, we obtain that for a typical configuration, the probability of k nodes being present in the Csr is 

_ ( c iV™\ f c 3 \ k A c 4 Y 2V ^~ k _ (cic 3 ) fc exp(-C2C 4 ) 



for large enough n. If Csr has k nodes not including S and R then S certainly has no more than k + 1 neighbors. 
In this situation, R is chosen by S with probability at least l/(k + 1). Further there are at most k other source 
nodes and the probability that no other node in Csr becomes active is at least (1 — pA) k - Thus, 

PmEl) > g (o 1 c^exp(-c 2 c 4 ) _|_ (1 _ pA)t 



k=0 



^ (C1C 3 (1 - PA)) 
> exp(-c 2 c 4 )g (fe + 1) , • 

It is easy to see that for < pa < 1, the term on the right hand side is a lower bounded by a strictly positive 
constant. Hence, P(E2\E\) is strictly positive. This completes the proof of the lemma. ■ 

Theorem 1: Scheme IT achieves T(n) = 0(1). 

Proof: Consider a typical configuration so that policy S n is used. Fix a source node S and a relay node R. 
Let A(t) be the number of bits transmitted from S to R in sub-slot A of time-slot t. If S transmits to R successfully 
in sub-slot A of time-slot t, A{f) = W/2 otherwise A(t) = 0. 

First we determine 2£ [>!(£)]. Let Fi be the event that S and R are neighbors and F 2 be the event that S transmits 
to R successfully. Then 

W W 
E[A(t)) = ~^P{Fi n F 2 } = —P{F 1 }P{F 2 \F l }. (3) 

From Lemma |3a), P{F 2 \Fi] > C5 > 0. Due to the independent motion of nodes S and R according to natural 
random walks, the joint description of their positions is a two-dimensional random walk on a discrete torus of 
size y/n x y/n. It is easy to see that the stationary distribution for this process is the uniform distribution on n 
joint positions. Since S and R become neighbors when they are in one particular joint position out of these n joint 
positions, it follows that the probability of S and R being neighbors is 1/n, i.e., P(F\) = 1/n. Hence from © it 
follows that E[A(t)] = 0(l/n). 
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Now the positions of nodes S and R form an irreducible, finite state Markov chain and A(t) is a bounded, 
non-negative function of the state of this Markov chain at time t. Therefore by the ergodicity of such a Markov 
chain, the long-term throughput between S and R is 

1 T 

k-^^^Hefi/n). 

t=i 

Thus the throughput between a source node S and any other node in sub-slot A is 8(l/n). Similarly, it can be 
shown that the throughput between any node and a destination node D in sub-slot B is also 0(l/n). The value of 
< a < 1 guarantees that the arrival rate of packets belonging to every S-D pair at any relay node is strictly less 
than the service rate. This ensures the stability of the queues formed at the relay nodes, which in turn implies that 
the throughput between each S-D pair is simply the sum of the throughputs between S and the other n — 1 nodes 
in sub-slot A. Hence the throughput of each S-D pair is 6(1). 

We have shown that in a typical configuration, Scheme IT provides 6(1) throughput between all S-D pairs. From 
Lemma ^ configurations are typical whp. Hence it follows that Scheme n has throughput T(n) = 6(1). ■ 

Note that for the unrestricted mobility models in [8] and [6], it is possible to prove a stronger result that each 
S-D pair has 6(1) throuhgput for any n with probability 1 instead of probability approaching 1 as n tends to 
infinity, as in the present case. 



IV. Delay of Scheme II 



Under Scheme IT, if the configuration is not typical, direct transmission is used, in which case the delay for each 
packet is 1. Since the delay of a scheme is defined to be the expectation over all configurations of the average 
delay, the delay for Scheme IT is determined by the expected delay over typical configurations. So we shall assume 
that the configuration is typical. 

Consider a particular S-D pair. Packets from S reach D either directly by a single hop in sub-slot A or through 
any of the other n — 2 nodes, which act as relays. Since the nodes perform independent random walks, only 6(l/n) 
of the packets belonging to any S-D pair reach their destination in a single hop. Thus, most of the packets reach 
their destination via a relay node, in which case the delay is two time-slots for two hops plus the mobile-delay, 
which is the time spent by the packet at the relay node. 

Each relay node maintains a separate queue for each of the S-D pairs. Fix a relay node, R, and consider the 
queue for the S-D pair under consideration. The mobile-delay mentioned above is the delay at this relay-queue.To 
compute the average delay for this relay-queue, we need to study the characteristics of its arrival and potential 
departure processes. 

First we obtain a lower bound on the delay at the relay-queue. Each node performs a random walk on a 1-D 
torus of size y/n on its great circle. We say that an S-D pair intersects node R's great circle k vertices apart if the 
lattice points where R can become neighbors of S and D are k lattice points (vertices) apart on the 1-D discrete 
torus of R. 

Fix an S-D pair and consider a particular relay node R. When a packet is transmitted successfully from S to R, 
D is equally likely to be in any of its yfn lattice points since it performs an independent random walk. Let Ty be 
the random time it takes for a random walk on a y/n x yfn torus to hit (0, 0) starting from (i, j). If the S-D pair 
intersects the great circle of R i vertices apart then the expected delay for packets of this S-D pair relayed through 

R is lower bounded by Ylf^Q 1 ^ij- 

Using the Chernoff bound for the sum of i.i.d. Bernoulli random variable (e.g., see [11]), it can be shown that 
Q(y/n) S-D pairs intersect the great circle of each node i points apart for < i < yfn — 1 whp. Hence the delay 
of Scheme IT, which is the expected delay over all packets is 



D(n) = O \E 



n ' 



As shown in [1], E 



Ti 



0(nlogn). Therefore, 



D(n) = fi(nlogn). 



(4) 
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The rest of this section derives an upper bound which is of the same order as the lower bound. It is hard to 
obtain an upper bound on the delay in the relay-queue since the arrival and service processes are complicated and 
dependent. We progressively obtain queues that are simpler to analyze and upper bound the delay of the previous 
queue as follows. We first upper bound the delay in the relay-queue by that in another queue, Q\, in which the 
arrival process is simpler. The delay of Q\ is upper bounded by that in Q2, which has a relatively simpler service 
process. However, the arrival and service process are not independent. The final part consists of introducing a virtual 
server with i.i.d. Geometric service times to break this dependence. With this overview, we proceed to the details. 

Recall that a packet arrives at the relay-queue when (i) S and R are neighbors, (ii) S becomes active (which 
happens with probability ap&), (iii) S chooses R as a random relay, and (iv) the transmission from S to R is 
successful. Similarly, a packet can depart from the queue when (i') R and D are neighbors, (ii') R becomes active 
(which happens with probability p&), (iii') R chooses D as the destination node, and (iv') the transmission is 
successful. We call such a time-slot a potential departure instant and the sequence of inter-potential-departure times 
is called the potential-departure process. Let the potential-departure process of the relay-queue be called {Si}. The 
qualifier potential is used since a departure can occur only if R has a packet for D. 

Consider a queue Q\ in which arrivals happen whenever (i), (ii) and (iii) above are satisfied, irrespective of 
whether (iv) is satisfied or not. The potential departure process for Q\ is the same as that for the relay-queue. Then 
it is clear that the expected delay in Q\ provides an upper bound on that in the relay-queue. 

Recall that the motion of each node is an independent 1-D random walk on a discrete torus of size yfn. We will 
say that two nodes meet when they become neighbors. Since nodes move independently the joint position of nodes 
R and D is a random walk on a y/n x y/n discrete torus and R and D become neighbors when the 2-D random 
walk is in state (0, 0), without loss of generality. Therefore, the inter-meeting time of R and D is distributed like the 
inter-visit time of state (0, 0) of a 2-D random walk. Since this is a Markov chain with n states having a uniform 
stationary distribution, we know that the sequence of inter-meeting times of nodes R and D, denoted by {n, i > 0}, 
is an i.i.d. process. Further, if r is a random variable with the common distribution then 

E[t] = n. (5) 

However a potential departure instant does not occur each time R and D meet. A potential departure instant 
occurs only if R also becomes active, chooses D as the random destination and the transmission is successful. If 
R and D are not chosen in spite of being in the same cell, it increases the likelihood of there being many more 
nodes in the same cell. Due to the random walk model of the node mobility, if there is a crowding of nodes in 
some part of the network then it remains crowded for some time in the future. Hence due to the Markovian nature 
of node mobility, the inter-potential-departure times are not independent. 

We want to obtain an upper bound on the delay of Q\ which has potential-departure process {Si}. To do this 
we will consider a queue, Q2, which has the same arrival process as Q\ but a different departure process {Si} 
such that Si < Si. Then the expected delay in Q2 would provide an upper bound on the the expected delay in the 
relay-queue. 

Nodes R and D perform independent random walks on 1-D tori of size y/n on their great circles as shown in 
Figure El and R and D meet when both are at a particular pair of lattice points. This is represented schematically in 
Figure |3j where R performs a vertical 1-D random walk and D performs a horizontal 1-D random walk. The joint 
motion of nodes R and D is equivalent to a random walk on a 2-D torus of size y/n x y/n and R and D meet when 
this 2-D random walk is in state (0, 0). The inter-meeting times of nodes R and D correspond to the i.i.d. process 
{7-;}. Further, let aj = n + . . . + t\ for % > 1, i.e., ai is the time-slot in which R and D meet for the ith time. In a 
typical configuration, we know that the number of other great circles that pass through Crd is © (y/n)- Allowing 
for the worst case, based on Lemma [0 let there be c\y/n = m — 2 other great circles that pass through Crd- These 
can also be thought of as performing independent random walks on the horizontal 1-D torus. Let nodes R and D 
be numbered 1 and 2 and the other c\y/n nodes be numbered from 3 to m and let X(t) = (X\(t), . . . ,X m (t)) 
denote the position of these m nodes on the y/n x y/n discrete torus at time t. 

A constant number of lattice points of the 1-D torus correspond to Crd and these are shown by the shaded region 
in Figure [3] and is referred to as set A. Let Ei be the indicator for the event that R chooses D and the transmission 
is successful in time-slot ctj. That is, Ei is the indicator for the event that a.\ is a potential departure instant. Let 
Ni be the number of other destination nodes in A in time-slot a». Then P{Ei = 1} depends on Ni only. Now, 
Ni depends on X(cti) which depends on the past given by E 1 " 1 = {Eq, . . . ,E^i} and t % = {tq, . . . , Tj}. Thus 
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Fig. 3. Schematic representation of the motion of nodes R and D on their respective great circles with £ 



the potential-departure process is generated by choosing some of the meeting instants of R and D according to a 
probability modulated by Ni, which is another independent process as shown in Figure |4] 
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Fig. 4. The 'x' marks correspond to the times when R and D meet each other. At some of these meeting instants the R-D transmission 
can be successful. Such points have been circled and correspond to Ei = 1. The inter-potential-service times are thus the sum of a few of 
the inter-meeting times of R and D. 



Above we described how the process {Si} can be generated using the processes {N} and {r^}, which in turn 
were obtained from {X(t)}, which corresponds to the independent random walks of all m nodes. Next we shall 
perturb the process {X(t)} to obtain {X(t)} and the corresponding {fj} and {Ni}. Let Z(t) be a 1-D horizontal 
random walk on a torus of size yfn. Let Xi(t) = Xi(t) + Z(t) be the position of node i, 1 < i < m, where 
the addition is modulo y/n. Then the inter-meeting times of any two nodes are the same as before since the 
position of each node is shifted horizontally by the same amount due to Z(t). As a result the processes r, and f, 
are identical. Under the modified setup, the lattice point at which R and D meet can be any element of the set 
B = {(i, 0) : < i < y/n — 1} instead of always being (0,0). Similarly, let Ni be the number of other destination 
nodes in the set A + Z(t). Then, {NA- is identical to {Ni}. Thus the process {Si} can also be generated (through 
{Ei}) using {N} and {fj} instead of {N} and {rj}. Therefore we shall use X{(t) as the position of node i at 
time t instead of Xi(t). Under this perturbed motion, R can be seen as if it performs a 2-D random walk on the 
yjn x ^fn torus while D and the other m — 2 nodes perform a 1-D random walk on a 1-D torus of size yjn which 
is subset B of the 2-D torus. Moreover, given X™(ai) = (X^(ai), . . . ,X m (ai)), P{Ei = 1} is independent of 
everything else. 

Lemma 3: There exists a constant (independent of n) cq > such that 

p{Ei = llr^" 1 ) > c 6 > 0. 
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Proof: The initial position of R, X\ (0) has a uniform distribution of the yjn x sfn torus. The initial positions 
of D and nodes 3 to m have independent uniform distributions on subset B = {(i,0) : < % < yjn — 1} of the 
^fn x y/n torus. As a result Xi(a\) = X 2 (ai) = / where / is a random variable with a uniform distribution over 
B. 

Let V = {Xz{ai), . . . , X m (a>i)) be the configuration of the m — 2 nodes other than R and D. Then the conditional 
probability of a potential departure given the past can be written as 



P{E l -l\r,E ) - ^ pM pjj=ri 



( P(r%P* 



> minP(P; = l|I/,r*,£ 
= minP(Pi = 1\V,t\E 



i-lx 



= minP(P i = (6) 

where the last equality holds because Ei is independent of everything else given V. 

Given a congfiguration V, the number of nodes in A + [i — 1,0) for i = 1, . . . , -y/ra torus can be found and 
this in turn determines the P {Ei = V\V). Hence, if Vi denotes the number of nodes other than R and D in the set 
A + (i — 1, 0) for i = 1, . . . , yjn then we can equivalently let the configuration be V = (Vi, . . . , V/^). 

Now consider a fixed configuration, V = v = (yi, . . . , v^), and let Z be a random variable which takes value 
i>i, 1 < i < with probability l/y'ra. Let ^4 consists of c 2 (some constant) elements. Then 

Recall that X^ctj) = /, where / is a random variable with uniform distribution on B. Further, from the description 
of Scheme IT, if there are destination nodes other than D in Crd then Ei = 1 if R chooses D out of all destination 
nodes that are its neighbors and the other vt nodes do not transmit. Since Crd contains all neighbors and more, 
the number of neighbors can be no more than Xi and hence for k = 1, ... , \/n, we obtain 

P(P, = 1\V = w,Xi(ai) = (k — 1,0)) > PA(l-PAr* + \ (8) 

v k + l 

Define a real valued function / : R — > E where /(x) = PA ^~^ — . It is easy to check that /(•) is a convex 
function. Hence, by Jensen's inequality, 

E[f(Z)] > f(E[Z]). (9) 

Using (0, (|8jl and Q, for any configuration V with corresponding v, we obtain 

v/n 

P(^ = l|T/ = «) = ^P(E i = l|F = t; J Xi(a i ) = (fc-l,0))P(Xi(ai) = (*!-l,0)|F = «) 

fc=i 

= -=y;P(^ = l|F = « J Xi(a i ) = (A!-l J 0)) 

> i ^ PA(i-PAr fc+i 

= > /(^[Z]) 

Combining (|6ll and dlOt completes the proof of the lemma. 



/(^)^>o. 
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Recall that the process {Si} is generated from {t{\ and {Ei\. Consider an i.i.d. Bernoulli process {E{\ with 
P{E\ = 1} = C6- Now we can construct a process {Si,i > 1} similar to the process {Si} using {t^} and {Ei} 
instead of {Ei}. Lemma [3] shows that the processes {Si} and {Si} are coupled such that Si > Si (the inequality 
corresponds to standard stochastic dominance). Now consider queue, Q%, with the same arrival process as Q\ but 
with potential-departure process {Si}. Depending on the value of cq, the value of a can be chosen so that the 
arrival rate is strictly smaller than the potential departure rate in Q2 so as to ensure stability. The distribution of 
Si is the same as t\ + . . . + tq, where G is an independent Geometric random variable with parameter c%. As a 
result, for any r G N, 

E[S r 1 ] = e(E[r r 1 ]). (11) 

In light of (fTTl . it is easy to see that the delay scaling of queue Q2 is the same as the delay scaling of a queue 
in which an arrival happens each time S and R meet with probability 0.5 and a potential departure occurs each 
time R and D meet. Since we are interested only in the delay scaling, henceforth we assume that in Q2, an arrival 
happens when S and R meet with probability 0.5 and a potential departure occurs whenever R and D meet. 

At this stage we have upper bounded the delay in the relay-queue by the delay in Q2. The inter-arrival times 
and the inter-potential departure times in Q2 are i.i.d. processes. However these two processes are not independent 
for the following simple reason: if the S-D pair intersects the great circle of R, k > vertices apart then R has to 
travel at least distance k on the discrete torus after an arrival for a potential departure to occur. 

Next, we will bound the delay in Q2 by the sum of the delays through two virtual queues, Q3 and Q4, in 
tandem. Both Q3 and Q4 will be shown to have delay of 0(n log n). This will imply that the delay of Q2 is 
0(n log n). Queues Q3 and Q4 are constructed as follows. The arrival process of Q3 is the same as that of Q2. 
The potential-departure process of Q3 is an i.i.d. Bernoulli process with parameter 2/3n (or potential departure 
rate j-). An arrival occurs at Q4 whenever there is a potential-departure at Q3. If Q3 is non-empty, then the arrival 
to Q4 is the head-of-line packet transferred from Q3 to Q4 or else a dummy packet is fed to Q4. Thus the arrival 
process at Q4 is the same as the potential-service process at Q3. By construction, the delay of a packet through this 
tandem of queues, Q3 and Q4, upper bounds the delay experienced by a packet through Q2. Now, from Lemmas |5] 
and [6] stated in the next section, the expected delay through Q3 and Q4 is 0(n log n). Thus the expected delay of 
the packets of each S-D pair relayed through each relay R in a typical configuration is O(nlogn). The delay of 
scheme is the expectation of the packet delay averaged over all S-D pairs and all relay nodes. Hence it follows 
that the delay of the scheme is O(nlogn). Combining this with we have proved the following. 

Theorem 2: The delay of Scheme IT is 0(nlogn). 

V. Remaining Proofs 

In this section, we prove Lemmas [5] and |5J which were used to prove that Scheme IT has delay of 0(n log n). 
Before proving these, we present Lemma |4] which will be useful for both these proofs. 

Recall that each node performs an independent random walk on a 1-D discrete torus of size ^fn on its great circle. 
Let Z be a random variable which is distributed as the inter-meeting time of two distinct nodes. The following 
lemma provides the first two moments of Z. 

Lemma 4: 

E[Z] = n, E[Z 2 } = 9(n 2 logn). 
Proof: As nodes perform independent random walks, the joint position of two nodes is a 2-D random walk 
on the y/n x y/n discrete torus. Thus the inter-meeting time of any two nodes is equivalent to the first return time 
to state (0, 0) for this random walk on a ^fn x yjn torus. Since we are interested only in determining the exact 
order of the moments, we will consider a simple random walk. 

Let X(t) = (Xi(t), X2(t)) G {0, . . . , ^fn — l} 2 be a simple random walk on the ^fn x ^fn torus. Then the first 
return time to state (0, 0) is 

T = inf{t > 1 : X{t) = (0, 0), X(Q) = (0, 0)}. 

Note that X(t) is a finite-state Markov chain with a uniform equilibrium on the n states. For any finite-state Markov 
chain, the expectation of the first return time to any state is the reciprocal of the equilibrium probability of the 
Markov chain being in that state. Hence, E[T] = n. 

Define, T = M{t > 1 : X(t) = (0,0)}. Observe that T differs from T in that T is conditioned on starting 
at X(0) = (0,0). Let Eu^yTf k ^ denote the expected time to hit state (k,l) for the first time starting from state 
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Let .EVfTo] denote the expectation of To given that X(0) is distributed according to the uniform stationary 
probability distribution n. Then, 

i,j=0 i,j=0 

^/ri—l*fn—l ^/ri—ly/n—1 

= E E -<ij) E (i,j) T (k,i) = Yl Yl 7r ( i »iM fe .o £? (t I i) r (*,o 

i,j=0 k,l=0 i,j=0 k,l=0 

= nlogn, (13) 

where (fT2l holds because £V ■ E^^T^^ = • E^^T^n for any < fc, Z < -y/n — 1 due to symmetry of states 
corresponding to cells on the torus. For the validity of d 1 3I >. see page 11 of Chapter 5 in [1]. 
Using Kac's formula (see Corollary 24 in Chapter 2 of [1]) and d!3l ». we obtain 

E[T 2 ] = ° J .. = 2n 2 logn + n. 

7T(0, 0) 

(14) 

Therefore, we obtain E\Z\ = n and E[Z 2 ] = 6(n 2 logn). ■ 
Lemma 5: Let D3 denote the delay of a packet through queue, Q3, as defined above. Then, 

E[D 3 ] = 0(n log n). 

Proof: An arrival occurs to Q3 when S and R meet with probability 0.5. Let {X} be the sequence of inter- 
arrival times to this queue. Then, Xj are i.i.d. with E[X{\ = 2E[Z] = 2n and E[Xf\ = @(E[Z 2 ]) = 9(n 2 logn) 
from Lemma0] The potential-departure process is an i.i.d. Bernoulli process with parameter 1/1. 5n. Let {1^} be 
the sequence of service times then Yj is a Geometric random variable with mean 1.5n. Hence E[Yi] = 1.5n and 
^[Y] 2 ] = 0(n 2 ). By construction, the service process is independent of the arrival process and hence Q3 is a 
GI/GI/1 FCFS queue. Then, by Kingman's upper bound [12] on the expected delay for a GI/GI/1 - FCFS queue, 
the expected delay of Q3 is upper bounded as 

Em _ o(SM±35*}) _ pi2±=>) . 0( „io g „). < 15) 

■ 

Lemma 6: Let D4 denote the delay of a packet through queue, Q4, as defined above. Then, 

E[D 4 ] = O(nlogn). 

Proof: Consider the service process of Q4, which is 1 at a potential departure instant and otherwise. This is 
a stationary, ergodic process since the inter-potential-departure times are i.i.d. with mean n. The Bernoulli arrival 
process to Q4 is independent of the service process with mean inter-arrival time l.hn. Since the arrival and service 
processes form a jointly stationary and ergodic process with mean service time strictly less than mean inter-arrival 
time, the queue has a stationary, ergodic distribution with finite expectation as shown by [9]. Thus Q4 is stable. 

Let Qt be the number of packets in the queue in time-slot t and let Qi be the number of packets in the queue at 
potential departure instant i. Thus the process {Qi} is obtained by sampling {Qt} at potential departure instants. 
Let Ai + i be the number of arrivals between potential departure instants i and i + Then the evolution of Qi is 
given by 

Qi+i = Qi - 1{q,>o} + A+i- (16) 

Comparing the evolution of the process {Qi} with that of {Qt} shows that {Qi} also has a stationary, ergodic 
distribution. Let Z be the inter-meeting time of any two nodes as defined in the beginning of this section. Then 
since the arrival process is Bernoulli and the inter-potential departure times are i.i.d. with common distribution that 
of Z, it is clear the {A{} is a stationary process. Let Q, Q and A be random variables with the common stationary 
marginals of {Qt}, {Qi} and {A{} respectively. Then taking expection in ( 11 6t under the stationary distribution, we 
obtain 

P(Q>0)=E[A}. (17) 
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The arrival proces is i.i.d. Bernoulli and hence conditioned on Z, the distribution of A is Binomial (Z, 2/3n). Since 
E[Z] = n from Lemma we obtain 

E[A\ = E[E[A\Z]\ = E -4-1 = 2/3. (18) 

L 1 .077- 

Squaring (fTol i. taking expecation, using the independence of Qi and Ai + ± and then rearranging terms, we obtain 

2(1 - E[A])E[Q] = P(Q > 0) + E[A 2 } - 2E[A]P(Q > 0). (19) 
Using (fT7l and (fT8l in the above, we obtain 

L ^ J 1(1 — E[A]) 2 V 9 J 

Recall that conditioned on Z the distribution of A is Binomial (Z, 2/ 3n) and hence 

£[A 2 ] = £?[S[A a ]^]] = ^ + ^(£[Z 2 ]-£[Z]) 

I " 9^) + 9^ 6(n2l0gn) = 0(1 ° gn) ' (21) 
where we used Lemma |4] As a result it follows from (120b that 

£[Q] = G(logn). (22) 

Next, we will bound E[Q] using E[Q]. To this end, consider a time-slot t and let the number of potential 
departures before time-slot t be I(t). Thus time-slot t is flanked by potential departures I(t) and I(t) + 1. Then 
Qt < + Ar(t)+l- Also using the fact that {Qt} is ergodic, with probability 1, we have 

T 



E[Q] = lim 



T— >oo i * — ' 

fc=l 

1 i(r)+i 

i — >oo 1 * — * 

3=1 

/(T) + 1 1 /(T)+1 

= lim V 7 — V (Q j Z j+l + A j+l Z 1+l ) 

T^oo T I(T) + 1 ^ v 3 3 J J ' 

j=i 

= ^(e[Q]E[Z\ + ^E[^ (24) 

= O(logn). (25) 

We used the fact that I(T)/T—>1/E[Z] by the elementary renewal theorem [12] in d23l and the independence of 
Qj and Zj + i in (E3Ji. Now using Little's formula, since the arrival rate is 2/3n, we conclude that 

E[D A ) = E[A]E[Q] = yO(logn) = 0(n log n). 

■ 

VI. Conclusion 

In this paper, we studied the maximal throughput scaling and the corresponding delay scaling in a random mobile 
network with restricted node mobility. In [3], it was shown that a particular mobility restriction does not affect the 
throughput scaling. In this paper, we showed that it does not affect delay scaling either. In particular, we show that 
delay scales as D(n) = G(nlogn) for a network of n nodes, which is the same as the delay scaling without any 
mobility restriction. This was understood to be a consequence of the fact that in spite of an apparent restriction, 
essentially the node mobility remaining unchanged in the sense that (i) each node meets every other node for 
0(l/n) fraction of the time with only 0(1) other neighboring nodes, and (ii) the inter-meeting time of nodes has 
mean of 8(n) and variance of 0(n 2 log n). 
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